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A new approach to generalised Casimir type of problems is derived within the 
context of renormalisable quantum field theory (QFT). We study the simplest case 
of a massive fluctuating boson field coupled to a time-independent background 
potential. We use analytic properties of scattering data to compute the relevant 
Green's functions at imaginary momenta, which in turn yields a simple and efficient 
method to compute (one-loop) vacuum energy densities in QFT. Renormalisation 
is easily performed in the perturbative sector by identifying low order Feynman 
diagrams with the first few Born approximation to the Green's function. Numerical 
examples illustrate the efficiency of our approach. 



1 Introduction 

In this talk, I report on the technical aspects of a novel quantum field theory 
approach to generalised Casimir problems. The work presented here was done 
in collaboration with N. Graham, R. L. Jaffe, V. Khemani, M. Scandurra, 
O. Schroder and H. Weigel; for a more complete exposition, see refP. 

The traditional approach to Casimir systems imposes boundary condi- 
tions on a quantum field ab initio^ 1 . In reality, however, Casimir forces arise 
from interactions between the quantum field and matter, and no such inter- 
action is strong enough to enforce a boundary condition on all frequencies of 
the fluctuating field. To study whether an idealised boundary condition limit 
exists in which physical quantities become independent of the regularisations 
involved, wc have recently proposed^ to embed the Casimir calculation in a 
quantum field theory (QFT). The external condition is then replaced by a 
renormalisable interaction with a smooth background field that imposes the 
boundary condition in a certain limit. QFT renormalisation is the only sen- 
sible way to discuss and eventually remove UV divergences. Any remaining 
infinities in our approach indicate that the quantities under consideration de- 
pend in detail on the physical UV cutoff provided by the material. In ref3 
we argued that the traditional approach is perfectly valid for some quantities 
(e.g. forces between rigid bodies) while it fails for others (e.g. the Casimir 
stress on isolated surfaces). 

Since the boundary condition limit in the QFT approach involves arbi- 
trarily strong and sharply peaked background fields, the standard tools for 
the calculation of effective energies (or functional determinants) cannot be 
employed: Perturbation theory fails because of the strong coupling, while a 
derivative expansion is ruled out by the high Fourier modes in the background 
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field. Thus, we had to develop general new methods to compute renormalized 
one-loop quantum energies and energy densities. These methods are impor- 
tant in their own right and have a much broader application than the proper 
Casimir problem studied in re for details the reader is referred to refs. . 
In the present talk, I will concentrate on the derivation of our approach from 
conventional quantum field theory, which is presented in section [21 Compu- 
tational techniques for Green's functions and renormalisation are discussed 
in section |3| Section 0] presents a numerical example to demonstrate the effi- 
ciency of our method. I conclude with a brief summary and outlook on future 
directions. 



2 The Method 



For simplicity, I will restrict myself to the case of a massive scalar boson <j> 
interacting with a time-independent scalar background field er(x) through the 
coupling £ int = — -|0 2 cr(x). This model is renormalisable (in n < 3 space 
dimensions) and can be used as a QFT implementation of Dirichlet b.c. on 
a surface S^- if we let the background a be strong and sharply concentrated 
on iS. In general, our method requires that the background is sufficiently 
short ranged so that a conventional scattering theory can be defined. More 
importantly, we require enough (spherical) symmetry to separate the scatter- 
ing problem into partial waves; we shall consequently assume that cr(x) only 
depends on r = |x|. The vacuum energy density in a spherical shell of radius 
r is then 
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where |0) is the interacting vacuum and the operator in the vev is just the 
symmetric energy density operator Too for our model. To further evaluate 
eq. l(T)l. we first perform a partial wave decomposition of the field operator, 
</>(£, x) = 4>i{t-, r ) y{t}{x) where y{i}{x) are the n-dimensional spherical 
harmonics and {£} refers to the set of all angular quantum numbers in n space 
dimensions. Finally, we make a Fock decomposition for the radial operators, 



dk 
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+ b.s. 



_ , (2) 

with ui — V m 2 + k 2 . The wave functions ipe(k,r) in eq. J5J) are a complete 
set of scattering and bound state b solutions to the fully interacting field equa- 



a It should be mentioned that the world-line formalisnUl provides an alternative way for an 
exact (numerical) computation of 1-loop quantum energies. It is even applicable in cases 
where there is not enough radial symmetry for our method to work. 

'The term "b.s." in eq. indicates a similar contribution from the bound states of energy 



which are required for the completness of the function system. 
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tions. They can be ortho-normalised such that the expansion coefficients ag (fc) 
obey the standard commutation relations [ag(k), a^,(k')] — S(k — k')5ggi, (and 
similarly for the b.s.) which allows the usual particle interpretation (e.g. aj(fc) 
creates a V^-mode acting on |£1)). 

Inserting our decomposition back into eq. , the result is best expressed 
in terms of the Green's function 

Gt(r y. k) = .l r dg ^^faO_y^(0fcf) (3) 
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which is meromorphic in the upper complex fc-plane with simple poles at the 
bound state momenta k — in.j (in the next section, I will sketch an efficient 
method to compute Gg). The Green's function is intimately related to the 
density of states. In fact, defining a local spectral density 

pe(k,r) = -ikGg{r,r;k) , lmfc>0 (4) 

in the upper fc-plane, we easily recover the local density of states for real fc, 

Repi(k,r) = Im {fc Gg(r, r; fc)} = ip^(k,r)ipe(k,r) . 

The energy density, eq. Q, can be expressed (up to a total derivative) as an 
integral of the single particle energies u>(k) over all real fc, weighted by the local 
density of states in each channel £ (plus a contribution from the bound states 
ujj). For practical purposes, however, it is much more convenient to rotate the 
integration contour to the upper complex fc-plane. We get three contributions: 
(1) The residues from the bound state poles on the imaginary axis cancel 
the explicit bound state contribution in Gg exactlypl (2) the discontinuity of 
the square root cut in the single particle scattering energies u> = V m 2 + k 2 
yields an integral along the imaginary axis, k = it, t € [771,00], and (3) the 
big semicircle at infinity gives no contribution if the integrand falls off fast 
enough for \k\ — > 00. To ensure (3), we have to improve the decay of the 
local spectral density eq. J1J at large |fc|. Since the Born series to the Green's 
function converges at large |fc|, it is sufficient to subtract a few low order Born 
approximations. We use the notation 

[pl(k,r)] N = pe(k,r) ~ pf\k,r) p{ N) {k,r) (5) 

for the N times Born subtracted density (and similarly for the Greens func- 
tion). Formally, the Born series is an expansion of Gg in powers of the in- 
teraction, i.e. the coupling strength A. When used in the expression for the 
energy density e(r), the Born terms should thus correspond to the usual per- 
turbative Feynman series. Notice that we need a cutoff at this point since the 
low order Feynman diagrams and Born approximations to Gg are precisely 
the UV divergent parts of the calculation. Using dimensional regularisation, 
the identification of Born and Feynman scries has been established rigorously 
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at least for the lowest orderJ^. The final result for the energy density in a 
spherical shell of radius r is then c 



4(t 2 - m 2 ) 



N 



N 

+ J2^l + ^CT(r). (6) 

i=l 

Here, is the multiplicity in the channel £, D r = <9 r (<9 r — (n — l)/r) is 
a total derivative operator and eqt denotes the counter terms necessary to 
renormalise the N low order Feynman diagrams e™|j. The perturbative renor- 
malisation of €fd is standard and briefly sketched in the next section, where I 
also present an efficient method to evaluate the local density pg and its Born 
series. It should be emphasised that eq. © yields a finite renormalised vac- 
uum energy density for any smooth background <r(r). Renormalisation group 
arguments prove that the final result is cutoff- and scheme- independent. 

To compute the quantum energy E q [o~] we can simply integrate eq. 10 
over all radii. Using the formulcP 

dr[ Pi (k,r)] =i^liiF t (k) (7) 

valid in the upper complex fc-plane, we can relate the density of states in 
fc-space to the log-derivative of the Jost function Fg(k). Rotating back to the 
real axis from above, this turns into the well-known expression 

\pt(k)]o=Pt(k)- P ? ) (*0 = ^§ (8) 

for the (change in the) density of states in terms of the phase shift Sg. Thus 
we have the quantum energy (= change in the zero-point energy due to the 
background a) 

I -Ira V 'I' 

where vt(t) = In Fi(it) is the log of the Jost function on the imaginary axis. 

The continuation of eq. © back to the real axis has been used exten- 
sively in the pasP. Though eq. Q is superior from a computational point 
of view, the formula on the real axis has the merit of a clear physical inter- 
pretation in terms of one-particle states and densities. In the next section, I 
will briefly discuss the methods necessary to turn eq. © and © into efficient 
computational tools. 



c Since the vacuum bubble diagram is not inserted back, eq. JHJ really represents the 
change in the vacuum energy due to the interaction. 
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3 Computational Techniques and Renormalisation 



The Feynman series for the energy density is obtained from the usual pertur- 
bative expansion 11 

(Q|Tbo(x)|0) = ^Tr [f x (-d 2 - m 2 - a)- 1 ] (10) 

where T x is the coordinate space operator corresponding to the insertion of 
the energy density JIJ at the space-time point x. It has pieces of order cr° and 
a 1 which read, in momentum space, 

(k'\f^\k) = e i{k '~ k > [fc ,0 fc°+k'-k + m 2 ] (11) 

{k'\f^\k) = a{x)e l(k '- k)x . 

The Feynman series consists of all graphs with a single </>-loop and arbitrary 
insertions of Tx°\ tJ 1 ^ and <j(x) (from expanding the propagator in eq. I|1Q(| S |. 
To order cr 1 , for instance, we have two diagrams: (a) a single insertion of Tj 1 ^ 
and (b) one insertion of a and one insertion of Tx°\ Only diagram (a) (the 
tadpole graph) is divergent, 



iTr 

2i 



d d k 1 



X f\-d 2 m 2 )- 1 = >(x) 7^ -r^-s . (12) 



2° v ~' J {2Ti) d fc 2 

It may easily be renormalised (and even cancelled completely) by a counter 
term proportional to er(x). For the present model the tadpole counter term 
ci<r and one Born subtraction N = 1 are in fact sufficient for space dimension 
up to n — 2. In n = 3, an additional mass counter term C20 2 is required 
and we have to perform N = 2 Born subtractions. The coefficient C2 may be 
fixed by a specific choice of renormalisation condition but the renormalisation 
group ensures that physical quantities are scheme independent. 

More generally, the counter terms are local monomials in a and da. It is 
then immediately clear that renormalisation will only affect the energy density 
on the support of the counter terms. For a Casimir type of calculation, the 
background becomes sharply concentrated on the b.c. surfaces and the energy 
density away from the surfaces cannot be affected by the counterterms; it 
must thus be finite. The divergence of the energy density as one approaches 
Casimir surfaces has in fact been known for a long time!^ 

The final ingredient in our main formula, eq. Ipjl. is the appropriate 
Green's function, which may be expressed as a product of Jost (ft) and regular 
(4>i) solution to the scattering problem, 

G,(r,r';fc) = (■^^ f <)y^, 

Fi{k) 



d For the total energy eq. J^J, the Feynman diagrams are more easily computed from the 
functional determinant representation E q = ^ lndet(— □ + m? + cr) in euclidean space (T 
is the euclidean time interval). 
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In order to evaluate this formula on the imaginary /c-axis, it is important to 
cancel all oscillating pieces in the numerator, which would otherwise become 
numerically intractable in the upper complex fc-plane. We cancel the free 
Jost solution wi(kr) (a modified Hankel function) explicitly by the ansatz 
ft — we ge and 4>e ~ he/we- The result is a representation of the form 6 

V ' ' ' (2£-2 + n)g e (k,0) V ' 

in terms of two new functions ge and hf . Since both the Jost and regular 
solution obey the full radial field equation with certain boundary conditions, 
it is easy to derive ODE's for the relevant functions ge and he- The explicit 
expressions are somewhat lengthy and the reader is referred to refP for de- 
tails. A few general remarks are in order. By introducing two functions with 
boundary conditions at r = (for he) and at r — oo (for ge) we can avoid 
expensive shooting methods for the Green's function Ge- The solutions ge 
and he can be shown to be smooth bounded functions of r in the upper com- 
plex fc-plane. The ODEs may, in fact, be solved directly on the imaginary 
axis k = it, where the solutions are real and easily amenable to numerical 
computations. The Born approximations are most conveniently computed by 
iterating the ODEs for ge and he according to an expansion of Ge in powers 
of the interaction a. 

Finally, a similar ansatz known as the variable phase approach can be 
used for the log of the Jost function required in the quantum energy formula, 
cq. J2J. Writing ge(it, r) = e^ t,r ^ with a real function fy, we find the relevant 
function ve(t) = (3(t, 0) by integrating the ODE 

- #'(t, r) - [$(*, r)} 2 + 2t U(tr) #(t, r) + a{r) =0 (14) 

inwards, starting at r — > oo with j3e{t,r) = (3' t {t,r) = (Born approximations 
follow again by iteration). The coefficient function ^ is just the log-derivative 
of the free Jost function evaluated on the imaginary axis; for details see refffl 



4 A Numerical Example 

As a numerical example for the efficiency of our method, I discuss a back- 
ground field o~(r) which takes the profile of a Gaussian ring in n = 2 space 

dimensions, a(r) — Aexp[— 2lu ° ] • ^ n ^ ne where the width w tends to 
zero (and the strength to infinity) this background can be used to study the 
Casimir stress on a Dirichlet circle in n = 2 dimensions. Here, I will em- 
ploy it as a numerical example of how our method handles sharply peaked 
backgrounds. 

The left panel of figure ^ shows the contributions of various angular mo- 
mentum channels to the quantum energy density. These terms correspond 

e The s-wave in n = 2, 3 requires a slightly different treatmenffl 
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Figure 1. Left panel: Contributions of various angular momentum channels to the quantum 
energy density of the Gaussian ring. The first order Feynman diagram depicted in the right 
panel was excluded on the left. All quantities are measured in units of the boson mass m. 



to the once Born subtracted i-integral in eq. ©, i.e. they do not include the 
1 st order Feynman diagram-' 7 depicted in the right panel of fig. From the 
scales at the axes, it is clear that the first order diagram dominates the energy 
density (i.e. perturbation theory is accurate) even at relatively large couplings 
(fig. corresponds to A/m = 3.0). 

Figure [3 shows the complete quantum energy for Gaussian circles of vari- 
ous widths w. As can be seen from the height of the central peak, the density 
on the surface diverges in the sharp limit w — ► even after renormalisation. 
Our method allows to study this limit numerically by going to very small 
widths for which traditional approximation schemes fail. 

5 Conclusion 

In this talk, I have presented a new approach to the computation of 1-loop 
vacuum energies and energy densities in quantum field theory. Starting from a 
conventional Fock decomposition or, alternatively, from a Green's function ap- 
proach, we identify the potential UV divergent contributions to the quantum 
energy with the low order Born approximations to the Green's function. We 
subtract these contributions and add them back as Feynman diagrams which 
allows for a conventional and scheme independent renormalisation using stan- 
dard counter terms. The resulting formulae for the Casimir energy (density) 
are finite functionals of smooth background fields and easily amenable to nu- 
merical treatment. As an example, I have shown results for the Casimir energy 
density of a Gaussian circle in n — 2 space dimensions. 

Our method has many interesting applications. Among the topics cur- 
rently under investigation are magnetic vortices (which play a role in the 

■^Note that the tadpole diagram has been cancelled completely by the counterterm, but 
there is still a second (finite) contribution of order a 1 as explained in the last section. 
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Figure 2. Quantum energy density for Gaussian rings of various widths w. 



confinement properties of gauge theories) and Z-strings in electroweak the- 
ory, which bind fermions efficiently and may thus be of relevance in scenarios 
of baryogenesis. 
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